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Abstract 

We address the question whether a graviton could have a small nonzero mass. The 
issue is subtle for two reasons: there is a discontinuity in the mass in the lowest 
tree-level approximation, and, moreover, the nonlinear four-dimensional theory of a 
massive graviton is not defined unambiguously. First, we reiterate the old argument 
^ ■ that for the vanishing graviton mass the lowest tree-level approximation breaks down 

since the higher order corrections are singular in the graviton mass. However, there 
can exist nonperturbative solutions which correspond to the summation of the singu- 
lar terms, and these solutions are continuous in the graviton mass. Furthermore, we 
study a completely nonlinear and generally covariant five-dimensional model which 
mimics the properties of the four-dimensional theory of massive gravity. We show 
that the exact solutions of the model are continuous in the mass, yet the perturba- 
tive expansion exhibits the discontinuity in the leading order and the singularities in 
higher orders as in the four-dimensional case. Based on exact cosmological solutions 
of the model we argue that the helicity-zero graviton state which is responsible for 
the perturbative discontinuity decouples from the matter in the limit of vanishing 
graviton mass in the full classical theory. 



1 Introduction 



Could a graviton be massive? The naive answer to this question seems to be positive. 
Indeed, if the graviton Compton wavelength, \ g = mZ 1 , is large enough, let us say 
of the present Hubble size, we should not be able to tell the massive graviton from 
a massless one. In fact, astrophysical bounds are even milder, X g > 10 24 cm |]] (see 
also Refs. [0). However, in general relativity (GR) the issue turns out to be more 
subtle. A dramatic observation has been made in Refs. [0-0] according to which 
predictions of massless GR, such as the light bending by the Sun and the precession 
of the Mercury perihelion, differ by numerical factors from the predictions of the 
theory with a massive graviton, no matter how small the graviton mass is. This 
discontinuity, if true, would unambiguously prove that graviton is strictly massless 
in Nature. 

The arguments of Refs. [0-0] were based on the lowest tree-level approximation 
to interactions between sources. In this approximation the discontinuity has a clear 
physical interpretation. Indeed, a massive graviton in four-dimensions has five phys- 
ical degrees of freedom (helicities ±2, ±1, 0) while the massless graviton has only 
two (helicities ±2). The exchange by the three extra degrees of freedom can be inter- 
preted in the limit m g — > as an additional contribution due to one massless vector 
particle with two degrees of freedom ( "graviphoton" with helicities ±1) plus one 
real scalar ( "graviscalar" with the helicity 0). The graviphotons do not contribute 
to the one-particle exchange — their derivative coupling to the conserved energy- 
momentum tensor vanishes. The graviscalar, on the other hand, is coupled to the 
trace of the energy-momentum tensor and its contribution is generically nonzero. 
It is what causes the discontinuity between the predictions of massless and massive 
theory in the lowest tree-level approximation. 

However, as was argued in Ref. |J, this discontinuity does not persists in the 
full classical theory. It was shown that the lowest tree-level approximation to the 
calculation of interactions between two sources breaks down when graviton mass is 
small. The next-to-leading terms in the corresponding expansion are huge since they 
are inversely proportional to powers of m g . Thus, the truncation of the perturbative 
series does not make much sense and all higher order terms in the solution of classical 
equations for the graviton field should be summed up. The summation leads to the 
nonperturbative solution which is continuous when m g — > 0. The perturbative 
discontinuity shows up only at large distances where higher order terms are small, 
these distances are growing when m g — > 0. In other words, the continuity is not 
perturbative and not uniform as a function of distance. 

A simple reason why one could expect the violation of the lowest tree-level ap- 
proximation is that it does not take into account the characteristic physical scale 
of the problem; while the nonperturbative calculation of the Schwarzschild solution 
does account for this effect. In the nonperturbative solution the coupling of the 
extra scalar mode to the matter is suppressed by the ratio of graviton mass to the 
physical scale of the problem. Hence, the predictions of the massive theory could be 
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made infinitely close to the predictions of the massless theory by taking small m g . 

The argumentation can be conveniently presented by considering the gravita- 
tional amplitude of scattering of a probe particle in the background gravitational 
field produced by a heavy static source. This amplitude has the following generic 
structure (note, that we use the flat metric r]^ u = diag(— 1, 1, 1, 1)): 



a(q 2 ) V»t'-b(<?) Wt 



l V 



hT{q) C oc " ' 2 - ^ e - , (1) 
F q l + m z g — le 

where = p^p v and t' = p'jp' v refer to the heavy particle with the four-momentum 
Pn = (M, 0) and to the light particle with the momentum p'^ correspondingly [], see 
Fig. 1. The form factors a(q 2 ) and b(q 2 ) are functions of the momentum transfer 
q 2 and are defined by two parameters: the graviton mass m g and the Schwarzschild 
radius vm — 2Gj^M of the heavy particle with the mass M. 




Figure 1: Scattering of the probe particle at the gravitational field of the heavy 
source. The bold circle accounts for summation of the higher order iterations over 
the nonlinearities in the classical equations. 

In the lowest tree-level approximation of the massive theory the form factors a 
and b are just constants and the unitarity (sum over five helicities) fixes their ratio, 
a = 3b, while the same unitarity with two graviton states (helicities ±2) in the 
massless theory gives a = 2b. Therefore, the discontinuity H-^J appears. However, 
this is only valid for the small momenta q m g (m g r^) -1 ^, for which the higher 
order corrections are small 0. In the coordinate space it means that the linear 
approximation becomes valid only at the distance 

{m g r M f 15 {2G N Mm g f"> 

r^r m , r m = — y - = y - , (2) 

m g m g 

which for the Sun is bigger that the solar system size (see discussions in the next 
section). 

lr ro avoid the confusion note that we use only as a kinematical structure of the vertices not 
implying that it is the energy-momentum tensor. 
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On the other hand, at q ^> m g {m g r M )~ 1 ^ , i.e., at shorter distances, r «C r m , 
we expect that the summation of higher orders || returns the relation a = 2b of 
the massless theory. In other words, nonperturbative summation should lead to the 
decoupling of the graviscalar from the heavy source for distances r<r m . 

What was not verified in Ref. [BJ is a matching of the nonperturbative solution 
at r r m with the exponentially decreasing linear solution at r ^> r m . It might 
happen indeed that the solution matches an exponentially increasing function in- 
stead Boulware and Deser in expressed their doubts about an existence of the 
large distance matching. Moreover, they argued that there is no consistent inter- 
acting theory of the massive spin-2 field in 3+1 dimensions. One of the arguments 
in Ref. was that at quantum level the theory contained the sixth polarization 
in addition to the standard five polarizations. Furthermore, the mass term in the 
action is not uniquely defined beyond the quadratic order in the fields. 

These legitimate concerns can be addressed by embedding the 4D theory of a 
massless graviton into a five- dimensional theory — a route we take in the present 
paper. Indeed, gravity in five dimension is well defined as a classical gauge theory, 
a massless graviton has exactly five states. For the matter fields which are confined 
to the four-dimensional brane the theory mimics the massive spin 2 particle with 
the fifth component of the momentum playing the role of the massQ. 

The model which we discuss is that of Ref. ||. In this model matter is localized on 
a brane. The brane world-volume theory contains the induced 4D Einstein-Hilbert 
term due to which a five-dimensional graviton mimics the massive four- dimensional 
spin-2 state on the brane. In contrast with the 4-dimensional massive theory, in this 
case the full nonlinear action can be written. The two-body problem for sources on 
the brane is now well-defined. The amplitude has the same generic form ([!]) with 
substitution of q 2 + m 2 by q 2 + m c q, where m c is a counterpart of m g in the model. 
We present the arguments in favor of aforementioned behavior of the form factors 
a(q 2 ) and b(q 2 ). However, we did not manage to obtain the exact solution of the 
Schwarzschild problem in this case either. 

Instead, we derive a number of evidences supporting the conjectured behavior 
from the exact cosmological solutions P,|10| of the model. We show that the lowest 
tree-level perturbative result is off by a factor 4/3 as compared with the exact result 
and explain why the corresponding perturbation theory breaks down. Based on this, 
we expect that the perturbative discontinuity is absent on the nonperturbative level 
in the full classical theory indeed. 

Recently the problem of the vanishing graviton mass was studied in a different 



setup. It was shown in Refs. [11] and [12] that there is no mass discontinuity even in 



2 Such solution can still be acceptable as long as the exponential growth of the solution takes 
over at distances much larger than the observable size of the Universe. This will take place if the 
graviton Compton wavelength A g ^> 10 28 cm. 

3 Note the analogy with the supersymmetric BPS states whose mass is given by a central charge. 
This charge also can be viewed as an extra component of the momentum in the dimensionally 
enlarged space. 
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the lowest tree-level exchange on de Sitter (dS) |Tl],0 or Anti de Sitter (AdS) [ |TTi|T2" 
backgrounds []. This fits well with the discussions presented above. Indeed, in the 
case of the (A)dS background, even the lowest tree- level approximation does take 
into account the presence of a mass scale of the problem, which in that case is given 
by the cosmological constant A. It was shown in [EJ that the coupling of graviscalar 
is proportional to m 2 g /K when m g — > 0, and deviations from the massless model 
vanish in this limit. Since the cosmological constant in our world is restricted as 
A < 10 -84 GeV 2 , then the allowed graviton mass is in the range m g <C 10~ 42 GeV - 
i.e., the graviton Compton wavelength is bigger than our horizon size. The existence 
of such a tiny graviton mass is immaterial for all astrophysical and cosmological 



observations |TT| (see also an interesting discussion of the continuity issue in the 
recent work |TJJ ). Note, that the nonperturbative continuity allows for much wider 
range for the graviton mass, m g <C (7"A//r 5 ) 1//4 . Here r is the maximal distance from 
the Sun where the data are obtained, see Sec. 2 for numerics. 



In Ref. [15] it was argued that in the (A)dS background the perturbative discon- 
tinuity reappears at the one-loop quantum level — the phenomenon very similar to 
the one-loop discontinuity for a massive non-Abelian vector fields discussed in [|J. 
This is certainly true since the loops are sensitive to the number of particles running 
in the loop diagrams. From the practical point of view, however, the comparison 
of the theory with the experimental data on the light bending by the Sun and the 
precession of the Mercury perihelion is not affected by the small quantum loop cor- 
rections. Indeed, while the graviscalar decouples from the classical source it is still 
coupled to the graviton and does contribute to the quantum loops. However, such 
effects of quantum gravity are suppressed and most likely cannot be disentangled in 
solar system measurements. For these reasons, in what follows we are focusing on 
the (dis) continuity in the classical theory only. 

The paper is organized as follows. In Sec. 2 we recall the essence of the graviton 
mass discontinuity found in Refs. |3|-|5| and discuss the results of Ref. || where it was 
shown that there is in fact the continuity in the graviton mass in the full classical 
theory. In Sec. 3 we introduce the five-dimensional nonlinear model which mimics 
the properties of a four-dimensional massive gravitational theory. We show that the 
perturbative discontinuity which is present in the lowest tree-level approximation 
disappears in the exact solution of the model. In Sec. 4 we discuss another exact 
solution of the nonlinear model which interpolates between the four-dimensional and 
five-dimensional regimes. We conclude in Sec. 5. 



4 The consideration for the dS space is a bit subtle since for m 2 g < 2 A/3 (A being the cosmological 
constant) unitarity is violated in the theory fl"3j |. 
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2 Preliminaries: Massive Graviton in 4D 



We will consider the following action for a massive graviton on a flat 4D background: 



2 

ml 



S m = M&J d*x^\g\ \ R+-f K - (^) 2 JJ , (3) 

where g^ u = T]^ + h^ u and the Planck mass M P1 is related to the Newton constant 
Gn as Mpj = 1/(167tGjv)- The mass term has the Pauli-Fierz form [|l(J], in quadratic 



in h^v terms it is the only form which does not introduce ghosts [17] . We imply that 
indices in the mass term are raised and lowered by the tensor rj^. If it were g^ v 
instead the difference would appear only in the cubic and higher in h^ v terms which 
are not fixed anyway; higher powers of h^ v could be arbitrarily added to the mass 
term. 

In order to see the presence of the discontinuity in the lowest tree-level approxi- 
mation let us compare free graviton propagators in the massless and massive theory. 
For the massless graviton we find: 



D %;ap{q) = Q V^aVu/3 + ^ V^Vua ~ \ VlwvS) ' ( 4 ) 



where only the momentum independent parts of the tensor structure is kept. By a 
gauge choice the momentum dependent structures can be taken to be zero. On the 
other hand, there is no gauge freedom for the massive gravity given by the action 
(§), and the propagator takes the following form: 

D^afiiQ) = (o VrnVvf) + \ fjuPVua - \ fj^Vap) ~ 7 , ^ 2 ~ , (5) 



where 



q 2 + m z — xe 



Vnu = Vnv + 2" • 6 



Note the l/m^, l/m 2 g singularities of the propagator. 

The difference in the numerical coefficients for the r\ lLV , r] a ^ structure in the mass- 
less and massive propagators (1/2 versus 1/3) is what leads to the perturbative 
discontinuity [|3|-|5|. No matter how small the graviton mass is, the predictions are 
substantially different in the two cases. The structure ([5]) gives rise to contradictions 
with observations. 

To see how this comes about let us calculate the amplitude of the lowest tree- 
level exchange by a graviton between two sources with energy- momentum tensors 
T^ v and T' a g (the sign tilde denotes the quantities which are Fourier transformed to 
momentum space): 



Ao = -Stt G N f, u DT* TLp = (fr - \ ffj T >ttv . (7) 
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In the massive case this amplitude takes the form: 



An = -Stt G N % v DZ' aP f' ap = [% v - - v,., fjj) f"» . (8) 



q 2 + m 2 g \ ^ 3 

In the relativistic normalization we are using T^ u = (p|T M ^|p) = 2p^p u at zero 
momentum transfer, q — 0. Suppose we take two probe massive static sources with 
masses Mi and Mi. Then only Too, Tg are non- vanishing and the lowest tree- level 
graviton exchange determines the Newtonian interaction, 

i; {r) = / d3 g ^ <** _ G N M,M 2 



(2tt) 3 4MiM 2 r 

V (r) - f^S-e*^- ~ -- GnMiM2 c— (9) 
mK ) ~ J (2vr)3 4M X M 2 ~ 3 r 1 ' 

Expressions (|^) and (0) give different results for the Newtonian attraction even in the 
range r < A 9 where one can neglect the exponential decrease. This difference can 
be eliminated by redefining the Newton coupling for the massive theory as follows: 

4 

Gn = - Gn , (10) 

where Gn is the Newton constant of the massless theory. For nonrelativistic prob- 
lems the predictions of the massive theory with the coupling rescaled by a factor 
3/4 at m g — > are identical to those of the massless theory with the coupling Gn- 

However, this is not enough to warrant the viability of the massive model. The 
relativistic predictions in the two cases are different For instance, the predic- 

tions for the light bending by the Sun are in conflict. At the classical level the trace 
of the energy-momentum tensor for light is zero. Therefore, the second term on the 
right hand side of Eqs. (0) and @ is not operative for light. Hence, the amplitudes 
^4 an d Am are identical in this case. However, we have established above that 
calculations in the massive theory should be performed with the rescaled Newton 
constant. Taking into account this fact, the prediction for the light bending in the 
massive theory is off by 25% 

We could certainly take an opposite point of view. Namely, do not rescale the 
Newton constant of the massive theory. In this case the predictions for the light 
bending in the massive and massless models are identical. However, the Newton 
force between static sources would differ by a factor of 4/3. 

The above considerations are based on the lowest perturbative approximation. 
The question is whether these results hold in the full classical theory. Normally, one 
would expect that for the solar system distances the lowest approximation is well 
justified. However, it was argued in Ref. |J that the approximation breaks down in 
the massive theory for relatively short distances. Since this breaking manifests itself 
in a rather interesting way we will briefly summarize the results of Ref. below. 
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To see the inconsistency of the perturbative expansion in Gn let us look (fol- 
lowing ||) at the Schwarzschild solution of (|j). We parametrize the interval for a 
massive spherically symmetric body as follows: 



ds 2 = -e v{p) dt 2 + e^ p) dp 2 + e^ p) p 2 {d9 2 + sin 2 fl , 



11) 



In the massless theory the function p is redundant due to the reparametrization 
invariance of the theory; it can be put equal to zero. However, in the massive case 
this gauge symmetry is broken and p is nonzero. Therefore, in order to compare the 
results in the massive and massless case one has to do the substitution: 



/> <>X P \ g 



exp(A) = y + YJ^j exp(a-p). 



(12) 



The standard Schwarzschild solution of the massless theory takes the following 
form: 



^Schw . r 



(r) = -A Schw (r) 



In 1 



Tm \ r M _ 1 f r M 
r J r 2 V r 



+ . . . 



p Schw (r) 



. 



(13) 



Here Tm = 2G^M is the gravitational radius of the source of mass M. 

Let us compare this with the perturbative in Gn solution of the massive theory 
obtained in Ref. ||. In the leading plus next-to- leading approximation in Gat the 
solution reads: 



v 



r 

1 r JL 

2 r 

1 r M 



1 + 



7 r M 



32 rrigr 5 
21 r M 



mj. r 5 



2 m 2 r 3 



1 + 



21 r M 



4 rrig r 5 



(14) 



We note the following peculiarities of results (14): 



In the leading order there is the finite discontinuity in the expression for A: 
the result of the massless theory in (|i~3"D differs from the result of the massive 
model by a factor 1/2. This is precisely the discontinuity which is seen in the 
lowest approximation. 

The next-to-leading corrections in (|TJj) are governed by the ratio rj\,//m^r 5 and 
are singular in m g . 

For any given distance r there is a value of m g below which the perturbative 
expansion in breaks down. 
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These results are in correspondence with the perturbative series for the scatter- 
ing amplitude described by Feynman graphs in Fig. 1. The leading terms in the 
expansions ( |T3D are given by the diagram of the first order in the source, i.e., the 
diagram with one cross. The singular in m g terms in the propagator (H) do not 
contribute in this order. In the next order (the diagram with two crosses in Fig. 1) 
we have two extra propagators which could provide a singularity in m g up to 1/m®. 
Two leading terms l/m % g and l/m|j do not contribute again so the result contains 
only the l/m^ singularity as in Eq. (|I4]). 

To demonstrate how badly the expansion in powers of Gn breaks down let us 
take the largest allowed value for the graviton mass, m g = (10 25 cm) -1 0] and 
calculate the correction to the leading result in the gravitational field of the Sun. 
We will find that at distances of order of the solar system size, i.e., at r ~ 10 15 cm, 
the next-to-leading corrections in (|14|) are about 10 32 times bigger than the leading 
terms. Therefore, this expansion is unacceptable. 

For a light enough graviton, however, a consistent perturbative expansion could 
be organized in powers of m g . In this case one finds ||: 

v(r) = - T -f + O (m 2 g sf^j , A(r) = T -f + O (m 2 g yfi^Pj , 

,(r) = ^+0(my), (15) 

where only the leading terms in Tm/t are retained. These expressions are valid in 
the following interval: 

r M « r « r m , r M = 2G N M , r m = . ( i6 ) 

m g 

For the gravitational field of the Sun this would correspond to the interval: 

3 • 10 5 cm < r < 10 21 cm, (17) 

where the lower bound is less than the radius of the Sun and the upper bound is of 
the order of a galaxy scale. Thus, for practical calculations within the solar system 
this expansion is well suited. 

As we see, the expressions for v and A in the leading approximation coincide with 
those of the massless theory (pT3|). Thus, there is no mass discontinuity. Moreover, 
the expressions ( |15|) explicitly shows non-analyticity in Gn, fJ> oc \/Gn, for v and A 
non- analytic terms are proportional to m 2 . 

We discussed in the Introduction subtle issues concerning the validity of the 
results discussed above arising even on the classical level: the nonlinear theory of 
massive gravity is not uniquely defined and it is complicated to make sure that the 
solutions which have no discontinuity do indeed satisfy the boundary conditions at 
infinity, i.e., that for r ^> l/m g the solution matches the exponentially decreasing 
function. 
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As we already noted even the exponentially growing solution can be acceptable 
when the graviton Compton wavelength becomes larger than the observable size of 
the Universe. The Yukawa factors due to the graviton mass, exp(±m 9 r), can be 
made to be arbitrarily close to the unity by decreasing the graviton mass. However, 
as we discussed above, this does not warrant the continuity of the m g — > limit since 
the coefficients in front of the perturbative potentials in the massive and massless 
theory (H) are different and m g independent. Therefore, the question whether the 
graviton could have a nonzero mass, effectively reduces to the question whether 
the graviton could have five polarizations. Indeed, these extra polarizations are 
responsible for the m g independent discontinuity in the coefficients in the potentials 
(§). Therefore, in what follows below we will address the question: "Can the graviton 
which describes the data in our observable Universe have five degrees of freedom?" 

In the next section we present a model based on five dimensions where massless 
graviton naturally has five degrees of freedom. The model is free of all the problems 
of the 4D massive gravity discussed above. We perform our analysis within this 
completely nonlinear theory in which exact solutions can be found. These solutions 
are compared with the perturbative results. We find that the picture outlined in 
the work || (and discussed above) holds. 



3 A Brane Model of Massive Graviton 

The 5D model which we will discuss was introduced in ||. The gravitational part 
of the action takes the form: 

S = Ml J d 4 x dy yJ\G\K+ J d A x ^\g~\ R(x) . (18) 

Where M* is a parameter of the theory and Mp\ = 1.7 • 10 18 GeV ^> M*. Fur- 
thermore, Gab is 5D metric tensor, A = {/i, 5} = {0, 1,2,3,5}, and 1Z is the five- 
dimensional Ricci scalar, g^ v denotes the induced metric on the brane which we take 
as 



9nv{x) = G^(x,?/ = 0), fi,v = 0,1,2,3 , (19) 

neglecting the brane fluctuations. 

We assume that our observable 4D world (4D matter) is confined to a tensionless 
brane (a tensionless hyper-plane in this case) which is fixed at the point y = in 
extra fifth dimension 0. In other words, we assume that the energy-momentum 
tensor of 4D matter has the following factorized form T^x) S(y). We also imply 
the presence of the Gibbons-Hawking boundary term on the brane, this provides 
the correct Einstein equations in the bulk. These simplifications help to keep the 
presentation clear and do not affect our main results. The brane world aspects of 
the model (|18|) were studied in detail in Refs. p|,|l8H20 . 



3 A simplest possibility is to consider a brane at a fixed point of the R/Z2 orbifold. 
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Let us study the gravitational potential between two static bodies located on 
the brane. This can be calculated from the action (|18|). The corresponding Green 
function is conveniently represented by working in momentum space in the four 
world- volume directions and in position space with respect to the transverse coordi- 
nate y. For the time being we can neglect the tensorial structure of the propagator 
(to be discussed below) and calculate the scalar part of the Green function. This 
can be done by calculating the corresponding propagator in a theory with scalars 
only which have the bulk and brane kinetic terms similar to . The result of the 
calculation reads as follows [fl: 



G(q,y = 0) 



where we introduce the parameter 



Mpj g 2 + m c ' 



m c = — 

r c 



2 Ml 



(20) 



(21) 



The Green function ([H]) has unusual features. It has a tachyonic pole at q 2 = 
~Qo + <f = m 1 which corresponds to the decay into the continuous tower of Kaluza- 
Klein states (which arise from the reduction of 5D graviton). Although, the five- 
dimensional graviton is well defined, from the 4D perspective it looks as unstable 
particle with the width m c . Nevertheless, the rules of integration for the propagator 
(p0|) in the complex energy plane can be defined consistently. 

In particular, using ( ^0|) we can find the static potential <fr(r). The result can be 
written in terms of special functions and has different asymptotic behavior for small 
and large distances (see Ref. ||). The "crossover scale" between these two regimes 
is defined by r c given in Eq. (pip. At short distances, i.e., when r<r c 



1 



1 f 7T 
- 77 + 



87r 2 M& r I 2 



•1+7 -In - 
r 



(22) 



Here 7 ~ 0.577 is the Euler constant. The leading term in this expression has 
the familiar 1/r scaling of the four-dimensional Newton law with a right numerical 
coefficient. The leading correction is given by the logarithmic repulsion term in (^D- 
Let us turn now to the large distance behavior. For r ^> r c one finds: 



16?r 2 M 3 r 2 



4 + of+ 



(23) 



The long distance potential scales as 1/r 2 in accordance with the 5D Newton law. 
Thus, the crossover scale (^T|) should be sufficiently large to avoid conflict with 
astronomical observations. In || it was estimated that for M* ~ 1 TeV, the crossover 
scale r c is around 10 15 cm, which is roughly the size of the solar system. This is 
too low to be consistent with data. Therefore, the scale M* should be taken to 
be at least a couple of orders smaller then 1 TeV. This is in no conflict with any 
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gravitational or Standard Model measurements (see discussions in Ref. |19],|2(J). We 
take r c > 10 25 cm which corresponds to M* < 1 GeV. 

The parameter m c plays a role in this model which in many respects is similar to 
that of the graviton mass m g in (|3|). Indeed, as m c — > 0, gravity on a brane becomes 
4D Newtonian at more and more larger distances. Moreover, the four-dimensional 
interaction in the model with the action ([Tj|) can be interpreted as an exchange of 
a four- dimensional state with the width equal to m c ||. In the next section we will 
find even more closer similarities between m c and m g . 



3.1 Perturbative Discontinuity 

To see that the model flT8"| ) exhibits the discontinuity in the one-graviton tree-level 
approximation let us calculate, following ||, the tensorial structure of one-graviton 
exchange. To this end we will solve the Einstein equations in the linear approxima- 
tion in Jiab which is the deviation from the flat 5D metric, 

Gab = Vab + h AB . (24) 
We choose the harmonic gauge in the bulk: 

d A h AB = \ d B h c c . (25) 

In this gauge from the {/i5} and {55} components of the sourceless equations of 
motion follows that 

/Vs = 0, h\ = h» . (26) 

Let us turn to the {/i^} components of the Einstein equations. After some 
simplifications they take the form: 

(Mld A d A + Mpi 8{y) d a d a ) V = - \t^ v - ^T a Q j 6{y) + 8(y) d»d v h a a . (27) 

There are two terms on the right hand side of this equation. The first one has a 
structure which is identical to that of a massive 4D graviton (or, equivalently of 
a massless 5D graviton). The second term on the right hand side which contains 
derivatives d fJi d u is not important at the moment since it vanishes when it is con- 
tracted with the conserved energy-momentum tensor. As a result, the amplitude of 
interaction of two test sources takes the form: 

KM.V = 0) f'^(q) oc ^ I " " , (28) 

q z + m c q 

where q = y/(p. We see that the tensor structure is the same as in the case of the 
massive 4D theory, see Eq. (H). 
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In analogy with the discussions in the previous section we could expect that the 
lowest tree-level approximation will break down in the next iterations in the classical 
source. Further indication on this is the existence of the singular in m c terms in the 
expression for the gravitational field h^ v produced by a static source. We write the 
energy-momentum tensor for the source as follows: 

T^(x) = -McV>W (3) (£), (29) 

where M is its rest mass. As before, let us make Fourier transform with respect to 
four world-volume coordinates. Then the solution looks as follows: 

h 00 (q,y) = c)-G N M——^- exp(-g|y|), (30) 

2 q 2 + m c q 

hij(q,v) = c\d N M 6tJ exp(-q\y\) 

4 q 2 + m c q 

+ cGN M^- r ^ exp(-q\y\), (31) 

m c q q l + m c q 

where c = — 167T. These expressions, taken at y — 0, should be contrasted with the 
lowest order expressions for the Schwarzschild solution in 4D theory with a massless 
graviton: 

h S oo hw (q) = c\G N M±- 2 , (32) 
z q 



h^(q) = c\g n M^. (33) 



Comparing the expressions (pO|-pl|) to those in (|32| ~ p3|) we draw the following 
conclusions: 

• Upon the substitution Gn — ► Gn the {00} components coincide for large 
momenta, or, equivalently for r ^ r c . 

• The {ij} component of the 5D theory consists of two terms. The first term, 
after the substitution G^ — > Gjy is twice as small as the corresponding term 
on the right hand side of the Schwarzschild solution fl33|). This is what gives 
rise to the discontinuity. 

• There is an additional term in the expression for hij(q,y = 0) which is pro- 
portional to: 

QiQj 
m c q ' 

This term does not contribute to the one-graviton exchange in the leading 
order because of conservation of the energy-momentum tensors (the diagram 
with a single cross in Fig. 1). However, it does contribute to higher order 
diagrams (the ones with two and more crosses in Fig. 1). This term is singular 
in m c and the perturbation theory in Gn breaks down when m c — ► 0. 
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Given these arguments, we conclude that for a consistent calculation of the 
interaction between two sources on a brane we should find the Schwarzschild solution 
which sums up all the orders of the Born expansion for the classical equations. 
Unfortunately, we could not manage to find the analytic solution. However, implying 
the existence of a smooth in m c — > limit, one could perform the expansion in m c 
in analogy with the 4D massive case |J . 

The {/iz^} component of the Einstein equation for the action ([18 ) can be inte- 



grated with respect to y in the interval —e<y<e with e — ► 0. The resulting 
equation takes the form: 



G lf u(x)+m c J G flv {x,y)dy = - - ^ S^SuoS^ (x) , (34) 



where Q^ u and Q^ u (x,y) denote the Einstein tensor of the worldvolume and bulk 
theories respectively. Since the extrinsic curvature has a finite jump across the 
brane, the second term on the left hand side of (|34T ) is nonzero even in the limit 
e — > 0. This term is proportional to the parameter m c with respect to which the 
expansion is performed (we imply that the metric is nonsingular in m c , this seems 
to be a reasonable requirement for a physically meaningful solution). 

Then, it is clear from (Rl) that in the lowest approximation in m c one recovers the 



usual 4D Schwarzschild solution of the massless theory fll3|) . For the calculation of 
the sub-dominant corrections in m c and for matching conditions at infinity, however, 
numerical simulations are needed. Note that in this case the solution should be 
matched at infinity to a well known 5D Schwarzschild solution which decreases as 
(r^f/r) 2 at infinity. This is an easier task compared to the 4D massive case where the 
power-low solution at short distances should be matched with the Yukawa potential 
at infinity 

Does this mean that we cannot compare analytically the perturbative and non- 
perturbative results in the model flIB"D? Not at all. Instead of finding the exact 
Schwarzschild solution we perform the similar analysis for other solutions which can 
be obtained explicitly. In the next section we discuss an exact nonperturbative cos- 
mological solution of the model found in Refs. |9|, |1(| which differs from the 
perturbative result by 4/3. 



3.2 Nonperturbative continuity 

In this section we study the cosmological solution in the model (|18f) found in Ref. || 
and ||10|| . It was already noticed in || that the cosmological evolution in (|i~8|) is 
governed by the Newton constant which differs from the the "Newton" constant of 
perturbation theory by 4/3. We will discuss in details this discrepancy. 

6 The next morning after this paper was submitted to the archive an interesting work pT] 
appeared. In Ref. pi]] the asymptotic form of the Schwarzschild solution for m c — > was also 
discussed and, moreover, certain generalizations of cosmological solutions of the model ( |l8| ) were 
obtained. 
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Our goal is as follows. We consider the solution of the model ( |T8| ) which describes 
the expansion of the matter dominated Universe. We will perform two distinct cal- 
culation for this. First we find the solution based on the Newtonian approximation. 
This calculation makes use of the lowest order potential between objects on the 
brane. As a second step we find the corresponding exact nonperturbative cosmo- 
logical solution of the Einstein equations. In the domain where the Newtonian 
approximation is legitimate, the perturbative result for the cosmological solution 
would coincide under the normal circumstances with the nonperturbative one as it 
happens in 4D world with a massless graviton. However, we find the discrepancy by 
a factor of 4/3 in these two methods. 

Let us start with the perturbative approach. As we established in the previous 
subsection the one-graviton exchange in the lowest approximation gives rise to the 
following expression for the potential of a massive source at short distances r -C r c : 

- M 

0(r) = -G N — . (35) 
r 

The appearance of the constant GV instead of Gn in this expression is related to 
the fact that we used the lowest tree-level approximation. 

Let us now use the standard consideration of the Newtonian cosmology[]. Con- 
sider a spherical ball with some uniform matter density in it. We assume that the 
radius of the ball R is much smaller then r c and that we are in a regime where 
the Newtonian approximation is valid. In this case the potential of the ball on its 
surface takes the form: 

M 

</>baii(r = R) = -G N — . (36) 

ri 

Let us consider a point-like probe particle of mass mo which is located just right on 
the surface of the ball. We neglect the back-reaction of this probe particle on the 
ball. The energy conservation condition for the system of the ball and probe takes 
the form: 

m R 2 ~ Mm km 

" G7V_ fT " (37) 

where dots denote the time derivative and k is some constant. We would like 
to calculate the time evolution of the radius R. In the regime which we discuss 
this is equivalent to the time evolution of the scale factor in Friedmann-Lemaitre- 
Robertson- Walker cosmology. In what follows we consider the solution which corre- 
sponds to the expansion of flat, i.e., k = 0, matter-dominated Universe. For k = 
we rewrite (1371) as follows: 



(38) 




7 For a careful treatment and interpretation of the Newtonian cosmology see, e.g., [E2| 
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where the density p for the matter-dominated Universe is related to the scale factor 
R as follows 



' = (39) 

where u is some constant. 

This is nothing but the Friedmann equation for the scale factor R for a flat 
matter-dominated Universe. We find the solution for the scale factor: 

R\t) = QiruG N t 2 . (40) 

This solution is consistent with the fact that we choose the time period when R ^ r c 
so that the brane world evolves in accordance with laws of 4D theory. What is 
important in our solution is the numerical coefficient in the relation ( |40"D which 
different from the 4D massless gravity case - it contains Gn = (4/3)Gat instead of 
Gn- Below we will show that the exact solution matches the massless gravity in the 
limit m c — > 0. 

Before discussing the exact solution let us explain why the Newtonian approach 
outlined above does not produce a correct coefficient. It is due to effects of nonlinear 
terms: similar to the Schwarzschild problem in 4D massive gravity discussed in 
Section 2 these corrections are defined by powers of the parameter 

Gu 1 

2 R 3 ~ ml t 2 ' { } 



It is clear that these corrections blow up at m c — > and we need to sum them up. 
The corrections seem to be small at the later time t ^> l/m c , but as we will see the 
4D approach stops to work at this epoch. 

Let us now solve the same problem using the exact Einstein equations. We 
parametrize the 5D interval in the following form: 

ds 2 = -N 2 (t,y)dt 2 + A 2 (t,y)dx i dx i + B 2 (t,y)dy 2 . (42) 

The 4D scale factor is defined as follows: 

R(t) = A(t,y = 0). (43) 



The solution was found in [91 and [ 1 C 1 : 



jV = l-|y||, A = R — \y\R, B = l, (44) 
R 

and the 4D scale factor obeys the following modified Friedmann equation: 

= y G N p - m c — . (45) 
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The m c — > limit of this equation is clearly incompatible with Eq. fl38|) which is 
based on the leading order approximation in the massive theory, but coincides with 
the result of massless gravity. This certainly implies that the Hubble parameter 
R/R is continuous in this limit — the assertion we verify below by presenting the 
exact solution of Eq. (|4"5j). 

We can absorb the parameters m c and Gn in Eq. (|45|) by rescaling, 




9 _ 



P 



3m c 
32nG 



N 



P 



■P, 
dp 

dr 



After introducing the variable 



x 



l + p= 1 + 



32nG 



N 



3 ml 



the exact solution can be written in terms of elementary functions for t(x) 



3 

- m c t 
2 



1 



1. yi+i 
,-T + 2 log 7i^T 



When r = m c i > 1 we get for the scale factor 

8tiGnu 



R* 



t 



m. 



1 _ log(3m c t) + 1 
3m r t 



(46) 



(47) 



(4J 



(49) 



This unusual (compare with the 4D Newtonian cosmology in Eq. fl4~0D) behavior is 
typical of a pure brane cosmology regime [23[] where one has if 2 oc p 2 - indeed, 
GN/m c = l/(327rM^) plays the role of GV m the 5D world. It is only relevant to 
the late time cosmology, t ^> l/m c — the epoch where the Hubble parameter is 
small, H ~ 1/t <C ?7i c , and the expansion enters the 5D regime, as analyzed in ||. 
Therefore, the 4D Newtonian cosmology is not applicable at this epoch. 



For t = m c t <C 1 



R 3 = QnG N ut 2 



3 

-m c t + 



(50) 



In correspondence with the discussed above difference of Eqs. (|45D and (pq) , we see 



that i? 3 at m c = is different from the expression in Eq. ([40]) which was obtained 
using the lowest tree-level approximation by the same factor 3/4 — it contains Gn 
instead of Gn- Note, that the exact expression for R 3 is linear in Gn — no higher 
orders are present. 

The exact solution considered above gives an explicit demonstration of the non- 
perturbative continuity in the limit m c — > 0. This continuity is not uniform — for 
the given value of t the parameter m c should be much smaller than 1/t. This is 
the most strong constraint on the graviton mass coming from cosmology, m c < H Q , 
where Hq ~ 10~ 42 GeV is the present day Hubble parameter. 
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4 An Interpolating Solution 



In this section we discuss a cosmological solution found in [|1(J and show that it 
interpolates between the regimes with 4D and the 5D tensor structures. 

Let us start with the brane action (PI) and in addition introduce in the theory a 
negative cosmological constant on the brane A b and the matter density p > \A b \ (we 
put pressure equal to zero for simplicity). The time evolution of such a 4D brane 
universe is interesting, it evolves asymptotically to a static Minkowski space on the 



brane without any fine tuning [IU|. The asymptotic form of the metric is as follows: 
ds 2 = -(1 + 6 \y\) 2 dt 2 + dx { dxi + dy 2 , (51) 
where the constant b is 

b = |A 6 |/4M, 3 . (52) 
In fact, this is a solution to the equation 

Kab -\g ab K= ^ T AB {x) 5{y) , (53) 

where the energy-momentum tensor on the brane is 

T M „ = diag (0, -A 6 , -At, -A 6 ) , T 5fl = T 55 = , (54) 

i.e. p + A b — > in this limit. To warrant the 4D behavior, the induced 4D Ricci 
scalar on the brane was added in |TT| . 

The important thing is that the early cosmology of this model is standard, with 
no discontinuity in the Newton constant. Indeed, the Friedmann equation is given 
in (|4"5|) where p should be substituted by p + A b . The Newton constant on the 
right hand side of this equation is the conventional 4D gravitational constant which 
reflects no discontinuity. This is true as long as the early cosmology is concerned. 

Let us now look at the late cosmology, more precisely at the form of the metric 
fl5"ip to which the solution asymptotes. The metric on the brane is Minkowskian and 
static everywhere with only dependence on y. For small values of y, which satisfy 
b\y\ <C 1, this metric can be obtained as a perturbation on the flat Minkowski space. 
Indeed, for small perturbations (|24]) in the harmonic gauge ( p5|) we find Eq. (|27|) 
with the energy momentum tensor defined in ([54]). This equation has the 5D tensor 
structure on the right hand side. Let us now notice that the energy-momentum 
tensor (|)4|) satisfies the relation: 

T ij ~ ^ Trjij = , i, j = 1, 2, 3 . (55) 

Therefore, the equation for hij is simplified. This is completely due to the 5D 
tensor structure; in fact had we have a 4D tensor structure, this would not be so. 
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Furthermore, the solution of equation ( p7|) in the gauge ( |25| ) can be written in the 
following form: 



'oo 



-h 



55 



[A 6 | 
2 M? 



1 2/1 







. 



(56) 



One can indeed verify that this solution coincides in the first order with the exact 
solution (|5lD . For this we perform the following gauge transformation of the exact 
solution (two different signs correspond to the two sides of the brane): 



y = sign(z) - 



1 + 2b\z\ + 2b 2 z 2 



1/2 



After this the metric takes the form 

ds 2 = -(1 + 26|;s| + 2b 2 z 2 ) dt 2 + dxUxi + 



[l + 2b\z\f 



l + 2b\z\ + 2b 2 z< 



dz 2 



(57) 



(5* 



which in the leading order coincides with the perturbative solution. 



Therefore, we conclude that the cosmological solution of Ref. |T0| does indeed 
provide an explicit example with both asymptotic regimes: at small distances (small 
Hubble radius) the behavior is 4-dimensional with the 4D tensor structure, whereas 
at large distances (large Hubble radius) the behavior has the 5D tensor struc- 
ture. In this sense the solution discussed above captures the important features 
of a Schwarzschild solution of 4D massive theory; this is not surprising since it is 
asymptotically (in time) Minkowski on the brane. 



5 Discussions and Conclusions 

We discussed a nonlinear five-dimensional generally covariant model which resembles 
many crucial properties of a massive graviton in four-dimensions. The mass disconti- 
nuity is present in the lowest tree-level approximation, however, this approximation 
breaks down for the vanishing graviton mass and all the tree-level graphs should 
be taken into account. The resulting expression of the nonperturbative classical 
calculation is continuous in the graviton mass. Thus, there is no mass discontinuity 
in the full classical theory. 

There are three extra degrees of freedom in the massive (or 5-dimensional) theory 
compared to the massless one. Among these degrees of freedom only the helicity 
state (the graviscalar) has a nonzero coupling to 4D matter. However, this coupling 
tends to zero in full classical theory as the graviton mass (or m c in the 5D example) 
vanishes. Thus, all the extra degrees of freedom decouple in the massless limit. 

The interesting issue which we did not discuss in the paper is the emission of a 
helicity gravitons. Based on our observations and using the unitarity arguments 
we expect that the nonperturbative amplitudes of the radiation of the helicity state 
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by 4D matter fields also vanish with the graviton mass, while they are non-vanishing 
in the lowest tree-level approximation as was shown in Ref. p4| . 



In the small mass limit the extra degrees of freedom of a massive theory form 
an independent sector which decouples from our matter as the graviton mass goes 
to zero. These degrees of freedom do interact with each other; moreover, in pertur- 
bation theory these interactions are singular in the limit m g — ► 0. Certainly, on top 
of the classical effects there is an issue of quantum loops which we did not discuss 
in the present work. However, the loop effects are suppressed and most likely they 
cannot be disentangled in existing measurements. 
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